Modeling urban street patterns 
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Urban street patterns form planar networks whose empirical properties cannot be accounted for 
by simple models such as regular grids or Voronoi tesselations. Striking statistical regularities 
across different cities have been recently empirically found, suggesting that a general and details- 
independent mechanism may be in action. We propose a simple model based on a local optimization 
process combined with ideas previously proposed in studies of leaf pattern formation. The statistical 
properties of this model are in good agreement with the observed empirical patterns. Our results thus 
suggests that in the absence of a global design strategy, the evolution of many different transportation 
networks indeed follow a simple universal mechanism. 

PACS numbers: 89.75.-k, 89.75.Kd, 89.65.Lm 
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Transportation networks -structures that convey en- 
ergy or matter from one point to another- appear in va- 
riety different fields, including city streets [ll, [2|, plant 
leaves^! , river networks [4] , mammalian circulatory sys- 
tems [sl , networks for commodities delivery [6] , and tech- 
nological networks [7] . The recent availability of massive 
data sets has opened the possibility for quantitative anal- 
ysis and modeling of these patterns and we focus here 
on the urban street network. Despite the peculiar geo- 
graphical, historical, social-economical mechanisms that 
have shaped distinct urban areas in different ways (see for 
example jsl and refs. therein), recent empirical studies 
0, i, M El, in, E, 113, [H have shown that, at least at 
a coarse grained level, unexpected quantitative similari- 
ties exist. The simplest description of the street network 
consists of a graph whose links represent roads, and ver- 
tices represent roads' intersections and end points. For 
these graphs, links intersect essentially only at vertices 
and are thus planar. Although the importance of net- 
works in geography and urban modeling has been recog- 
nized for a long time [16], comparably less attention has 
been devoted to generative models for planar graphs in 
the recent literature on complex networks [17]. Our aim 
is to propose a simple model for planar graph generation, 
based on plausible physical assumptions, and which re- 
produces several empirical findings. In the first part of 
this paper we discuss the empirical and quantitative sig- 
natures that characterize the topology of street patterns 
and which suggest the possibility to identify some gen- 
eral driving force steering the formation and evolution of 
street patterns. In the second part, we propose and dis- 
cuss a simple and parameter-free model based on a princi- 
ple of local optimality that quantitatively reproduces the 
above mentioned empirical features. The application of 
optimality principles to both natural and artificial trans- 
portation networks has a long tradition 0] and in most 
cases requires the minimization of a global cost function 
(such as the average total time for example), in sharp 



contrast with the model presented here. The rationale to 
invoke a local optimality principle in this context is that 
every new road is built to connect a new location to the 
existing road network in the most efficient way [l^ . The 
locality of the rule is implemented both in time and space 
during the evolution and formation of the street network, 
in order to reflect evolution histories that greatly exceeds 
the time- horizon of planners. The self-organized pattern 
of street emerges as a consequence of the interplay of the 
geometrical disorder and the local rules of optimality. 

In [l|, ^ measurements for different cities in the world 
are reported. Based on the data from these sources, we 
plot in Fig. [U the number of roads E (edges) versus the 
number of intersections N . The plot is consistent with a 
linear fit with slope ^ 1.44. When individual data points 
are considered, the quantity e = E/N = {k)/2 {{k) is 
the average degree of a node) shows a range of values 
1.05 < e < 1.69, in between the value e = 1 and e = 2 
that characterize tree-like structures and 2D regular lat- 
tices, respectively. These values are however not very 
indicative: planarity imposes severe constraints on the 
degree of a node and on its distribution which is gener- 
ally peaked around its average value. Few exact values 
and bounds are available for the average degree of clas- 
sical models of planar graphs. In general it is known 
that e < 3, while it has been recently shown [2^ that 
e > 13/7 for planar Erdos-Renyi graphs [20]. In Fig. [TJd, 
we plot the total length £ of the network versus N for the 
towns considered in py]. Data are well fitted by a power 
function of the form jj^N^ with /i ^ 1.51 and f3 ^ 0.49). 
The simplest hypothesis consistent with the data, at this 
stage, is that of an homogeneous and translational invari- 
ant structure. Indeed, a simple scaling argument, that 
could apply to a large family of planar graphs, includ- 
ing regular lattices, suggests that the typical distance 
£i between connected nodes scales as £i ~ where 

p = N/L'^ is the density of vertices and L the linear di- 
mension of the ambient space. This implies for the total 



FIG. 1: (a) Numbers of roads versus the number of nodes (ie. 
intersections and centers) for data from [1] (circles) and from 
(squares). In the inset, we show a zoom for a small number 
of nodes, (b) Total length versus the number of nodes. The 
line is a fit which predicts a growth as ^/N (data from 1]). 

length i ~ Eii ~ ^Ly/N. The discrepancies between 
the measured (k) and /i, given the error bars, are there- 
fore not enough to reject the hypothesis of an almost 
regular lattice. However, the network of roads naturally 
produce a set of non overlapping cells, encircled by the 
roads themselves and covering the embedding plane, and 
surprisingly, the distribution of the area A of such cells 
measured for the city of Dresden in Germany [15] has the 
form P{A) ~ with a 1.9. This is in sharp contrast 
with the simple picture of an almost regular lattice which 
would predict a distribution P{A) very peaked around il. 
The authors of [l5[ also measured the distribution of the 
form factor (j) = AA/ (ttD'^)^ (the ratio of the area of the 
cell to the area of the circumscribed circle) and found 
that most cells have a form factor between 0.3 and 0.6, 
suggesting a large variety of cell shapes, in contradiction 
with the assumption of an almost regular lattice. These 
facts thus call for a model radically different from simple 
models of regular or perturbed lattices. In the following, 
we describe a model where the set of 'centers' (repre- 
senting new homes, businesses, etc.) and the network of 
roads that connects them grow simultaneously. New cen- 
ters are introduced every rc > I time steps, and for the 
purpose of the present study we simply assume that the 
location of new points is given exogenously and we first 
assume them to be randomly and uniformly located over 
a square of given size. Finite segments (of fixed and small 
length) of roads are simultaneously added to the exist- 
ing network every r^^ = 1 < rc in order to account for 
the limited time horizon of planners. The algorithm that 
drives the construction of new portions of roads is based 
on a local optimality principle and aims at connecting to 
the network the still unconnected centers using as little 
as possible road-length. 

In order explain the algorithm, we illustrate it on the 
simple example of Fig.[2l We assume that at a given stage 
of the evolution, two centers A and B still need to be 
connected to the network. At any time step, each center 
can trigger the construction of a single new portion of 
road of fixed (small) length. In order to maximally reduce 
their distance to the network, both A and B would select 
the closest points Mi and M2 in the network as initial 
points of the new portions of roads to be built. If Mi and 
M2 are distinct, segments of roads are added along the 



FIG. 2: M is the closest network's point to both center A 
and B. The road will grow to point M' in order to maxi- 
mally reduce the cumulative distance A of A and B from the 
network. 

straight lines MiA and M2B. If Mi = M2 = M, it is not 
economically reasonable to build two different segments 
of roads and in this case only one single portion MM' of 
road is allowed. Our main assumption is that the best 
choice is to build it in order to maximize the reduction 
of the cumulative distance from the network (M) to A 
and B 

A = [d(M, A) + d(M, B)] - [d{M', A) + d(M^ B)] (1) 

The maximization of A is done under the constraint 
I MM' I = const. <C 1 and a simple calculation leads to 

MM' (XUA^UB (2) 

where ua and ub are the unit vectors from M in the di- 
rection of A and B respectively. The rule ([2j) can easily 
be extended to the situation where more than two cen- 
ters want to connect to the same point M. Already in 
this simple setting non trivial geometrical features ap- 
pear. In the example of Fig. [2] the road from M will 
develop a bended shape until it reaches the line AB and 
intersects it perpendicularly as it is commonly observed 
in most urban settings. At the intersection point, a sin- 
gularity occurs with ua ^ ub ^ ^ and one is then forced 
to grow two independent roads from the intersection to 
A and B. The above procedure is iterated until all cen- 
ters are connected. Interestingly, although the minimum 
expenditure principle was not used, the rule Eq. ([2]) was 
proposed by Runions et al. [3] in a study about leaf ve- 
nation patterns and we can follow their implementation. 
In particular, the growth scheme described so far leads 
to tree-like structures and we implement ideas proposed 
in in order to create networks with loops. Indeed, 
even if tree-like structures are on one side economical, 
on the other hand, they are hardly efficient (for example 
the path length along a minimum spanning tree scales as 
a power 5/4 of the Euclidean distance between the end- 
points [21]) and better accessibility is granted if loops 
are present. Following [3], we assume that a new cen- 
ter can trigger the construction of more than one portion 
of road per time step. An unconnected center s 'stim- 
ulates' the addition of new portion of roads from any 
vertex v of the road network (vertices correspond to any 
end points of the previously introduced road segments) 
that is in its relative neighborhood [22]. A node v be- 
longs to the relative neighborhood of s if for any node 
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FIG. 3: Snapshots of the network at different times of its 
evolution: for (a) t = 100, (b) t = 500, (c) t = 2000, (d) 
t — 4000 (the growth rate is here 77 = 0.1). At short times, 
we have almost a tree structure and loops appear for larger 
density values obtained at larger times (the number of loops 
then increases linearly with time). 



u (center or vertex of the road network) the inequality 
d{v,s) < mdix{d{s^u)^d{u^v)) holds [23], which captures 
the loosely defined requirement that v belongs to the rel- 
ative neighborhood of s if the region between s and v is 
empty. Centers can therefore be reached by more than 
one road, leading to the formation of loops. When more 
than one center stimulates the same point the prescrip- 
tion of Eq. ([2]) is applied and the evolution ends when 
the list of stimulated points is exhausted. 

The final road network is achieved starting from a 
small set of no centers connected by roads and iterat- 
ing the following two steps: (i) at every time multiple 
of Tc we add n new centers whose locations r are cho- 
sen randomly according to a given distribution P(r); (ii) 
the road network grows according to the algorithm de- 
scribed above. When a center is reached by all the roads 
it activated, it becomes 'inactive' and cannot stimulate 
the growth of a road any longer. We repeat (i-ii) un- 
til the total number of centers reaches a desired value. 
Although it is clear that the focus of the present pa- 
per is on the road network growth, it is important to 
stress that our model relies on a number of simplify- 
ing assumptions, the most relevant of which is the fact 
that the centers are independently located one from the 
other and from the structure of the road network. In 
fact, strong evidences [ijj, suggest that this is not 
the case, and integrating the correlation centers- network 
is the next most important step [25]. Despite this lim- 
itation, the model produces realistic results, in good 
agreement with empirical data (discussed below) which 
demonstrates that even in the absence of a well defined 
blueprint, non-trivial global properties emerge. In Fig. [3] 
example of patterns obtained for a spatially uniform dis- 
tribution of new centers are shown for different times. As 
time progresses, density increases, and the typical length 
from a center to the existing road network shortens and 
scales as £1 ~ l/y^ ^ I/aA as observed in the simu- 
lations. Beyond visual similarities with real cities, the 
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FIG. 4: Simulation results (averaged over 1000 configura- 
tions), (a) Total length of roads versus the number of nodes. 
The dotted Une is a square root fit. (b) Structure factor distri- 
bution showing a good agreement with the empirical results 
of TF] . (c-d) Rescaled distributions of the perimeter (c) and of 
the areas (d) of the cells displaying an exponential behavior. 



ratio e = E/N has initially a value around 1 (corre- 
sponding to a tree-like network) and increases very fast 
with N reaching a value around e 1.3 which is not far 
from the empirical finding (here and in the following, we 
checked that the results were robusts for different val- 
ues of the growth rate n/rc)- Consistently, in Fig [4^ 
we observe a relationship between the total length and 
N that is well approximated by a function of the form 
with a ~ 1.90, again in reasonable agreement with 
the empirical data. Panels b, c, and d of Fig. [4] show the 
collapse for different values of of the distributions of 
(/), A and the perimeter p of the cells, respectively. The 
excellent collapses show that the structures obtained are 
consistent with the hypothesis of homogeneity and trans- 
lational invar iance formulated above. We also note that 
the distribution of the (j) factor is peaked around 0.6 and 
essentially supported in the interval 0.4 < ^ < 0.7, in 
very good agreement with facts reported earlier [15]. A 
simple spatial uniform disorder and a plausible mecha- 
nism that connects the centers to the network can thus 
explain the non-trivial form of the factor distribution, 
but predicts an exponential behavior for the area distri- 
bution (Fig . [4ll), in disagreement with empirical obser- 
vations [l5|. In real cases however, the density of centers 
is not uniform. We therefore relax this assumption and 
assume, as supported by a previous empirical study [10], 
that the centers' distribution follows the population den- 
sity and decreases as P(r) = exp(— |r|/rc) where r is the 
distance from the central business district and I/tc the 
population density gradient. Although most quantities 
(such as {k) and £) are not sensitive to the centers' dis- 
tribution, the impact on the area distribution is drastic. 
Indeed, as shown in Fig. [5l we observe a power law de- 
cay with an exponent equal to 1.9 ± 0.05 in remarkable 
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agreement with the empirical result of [15] for the road 
network of the city of Dresden. This agreement confirms 

P(A) 




FIG. 5: (a) Network obtained for an exponential distribution 
of centers (1000 centers and rc = 0.1). (b) In this case, the 
area distribution is a power law (obtained for 5000 centers 
and 100 configurations). The solid line is a power law fit with 
an exponent 1.9 (for this size and with a smaller number 
of configurations, we observe fluctuations of this exponent of 
the order of 10%). 

that the simple local optimization is a good candidate for 
the main process driving the evolution of city street pat- 
terns but also shows that the center spatial distribution 
P(r) is crucial. 

More than 50% of the world population lives today in 



cities and this figure is bound to increase [2^. This mi- 
gration effect has dictated a fast and short-term planned 
urban growth which needs to be understood and mod- 
eled in terms of socio-geographical contingencies and of 
the general forces that drive the development of cities. 
Previous studies about urban morphology have mostly 
tried to identify specific mechanisms that have shaped 
distinct urban areas in different ways. Here we studied 
a simple model based on the assumption that road net- 
works develop trying to grant in an efficient and at the 
same time economic way connections to a set of 'centers'. 
The model accounts quantitatively for a list of descrip- 
tors that characterize the topology of street patterns, 
and in a more qualitative way for the tendency to have 
bended roads-even in absence of geographical obstacles- 
and perpendicular intersections. Interestingly, the opti- 
mality principle applied here turns out to be general and 
was implicitly at the basis of a model previously inves- 
tigated [3] to explain the formation of veins' patterns in 
leaves, pointing to an unexpected generality of the prin- 
ciple in the formation of transportation systems. This 
model is simple enough to allow many interesting gen- 
eralizations. In particular, our results suggest that the 
local optimality principle is a key ingredient for a more 
general model describing the co-evolution of the center 
distribution and the network [251] . 
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